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SUMMARY 

Buckling of a delaminated region can cause high interlaminar stresses 
which, in turn lead to delamination growth. Hence, buckling strain is an 
important parameter in assessing the potential for strength loss due to the 
delamination. The objective of this study was to predict the buckling of an 
elliptic delamination embedded near the surface of a thick quasi-isotropic 
laminate. The thickness of the delaminated ply group (the sublaminate) was 
assumed to be small compared to the total laminate thickness. Finite-element 
and Rayleigh-Ritz methods were used for the analyses. The Rayleigh-Ritz 
method was found . to be simple, inexpensive, and accurate, except for highly 
anisotropic delaminated regions. Effects of delamination shape and orienta- 
tion, material anisotropy, and layup on buckling strains were examined. 
Results showed that (1) the stress state around the delaminated region is 
biaxial, which may lead to buckling when the laminate is loaded in tension, 

(2) buckling strains for multi-directional fiber sublaminates generally are 
bounded by those for the 0 deg and 90 deg unidirectional sublaminates, and. 

(3) the direction of elongation of the sublaminate that has the lowest 
buckling strain correlates with the delamination growth direction. 


♦Research Associate Professor, Old Dominion University, Norfolk, VA. 
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INTRODUCTION 


Composite laminates often contain delarainations . Causes of delaraination 
are many and include tool drops, bird strikes, runway debris hits, and manu- 
facturing defects. The shape of a delamination generally resembles an ellipse 
[1,2]. Delaraination can significantly reduce the compressive strength and 
stiffness of the laminate, primarily because the delaminated region loses 
flexural stiffness. Even seemingly benign delarainations sometimes initiate 
localized buckling, thereby creating high interlaminar stresses and extensive 
delaraination growth. Delaraination growth may lead to structural instability. 

Buckling and postbuckling of a near surface delaraination has been 
addressed by several researchers both analytically and experimentally [1-7]. 
Analyses have been developed for strip [3,4], rectangular [1,5], circular [2] 
and elliptical [6,7] delaraination shapes. These configurations are shown in 
figure 1. As indicated in figure 2 the buckled region is referred to herein 
as the M sublarainate. ,, In refs. 1-7 , the laminate was always loaded parallel 
to one of the principal axes of the sublarainate. The buckling of an ellipti- 
cal sublaminate of arbitrary orientation and anisotropic material properties 
has not been addressed. 

The objective of this paper is, therefore, to investigate the buckling 
characteristics of a wide variety of elliptical anisotropic sublarainates 
delaminated from a quasi- isotropic base laminate. The thickness of the sub- 
laminate is assumed to be small compared to the base laminate. The sublami- 
nate is assumed to be made up of laminated fiber reinforced composite with 
fibers oriented at various angles relative to the loading axis. A general 
purpose finite-element program, STAGS [8], was employed to calculate the 
buckling strain of typical sublarainate configurations. Also, a simple energy 
method, based on a Rayleigh-Ritz (R-R) formulation, was developed for the 



parametric study* The R-R analysis was evaluated by comparing it to the 
finite element analysis. In this paper results are presented for variety of 
four ply sublaminates. 


NOMENCLATURE 
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half-length of an elliptic sublarainate, m 
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A 22 * A 26 » A 66 J 

b 

C 0 

C! >^2 

D 11 » D 12»°16 ( 
°22 * D 26 » D 66 j 

E 

G 

h 

IK] 

[K 1 MK 2 ],[K 3 ] 

[Kg] 

N 

N x » > Ny ? , N^tyi 

U 

u 

V 

v 

w 


x-y-z 


: f -y f -2 


inplane stiffness coefficients of the sublaminate, N/m 

half-width of an elliptic sublarainate, ra 
generalized displacement, m 
generalized displacements, ra ^ 

flexural stiffness coefficients of the sublaminate, N/ra 

Young 1 s modulus, MPa 

shear modulus, MPa 

thickness of the sublarainate, ra 

bending stiffness matrix of the sublarainate 

bending stiffness coefficient matrices 

geometric stiffness matrix of the sublarainate 

number of laminae in the sublaminate 

sublarainate stress resultants, N/ra 

strain energy of the deformed sublaminate, J 

displacement in the x f -direction, m 

potential energy due to loading, J 

displacement in y '-direction, m 

transverse (z) displacement, ra 

base laminate Cartesian coordinate system 

sublarainate Cartesian coordinate system 
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I 

i 


a 

c 


c 



0 

V lam 

n 


fiber angle of a lamina in a sublaminate, degree 

inplane strain 

buckling strain 

one-term buckling strain 

angle between x and x' axes 

Poisson's ratio of the laminate 

total potential energy of the sublaminate, J 


Subscripts 


i 

x »y 

x’,y’ 

Definitions 
Base laminate: 
Buckling strain: 
Delamination: 
Sublaminate: 


ply number 

along or transverse to fiber direction 
x or y direction 
x* or y' direction 

group of plies below the sublaminate 

base laminate strain at which the sublarainate buckles 
a separation between plies of a laminate 

a set of intact plies separated from a thick laminate by a 
delamination 


DESCRIPTION OF THE PROBLEM 

Figures 2(a) and 2(b) show plan and sectional views of a quasi-isotropic 
composite laminate with a single delamination. As shown, the x-y-z coordi- 
nate system is oriented such that the x-axis is along the laminate loading 
direction. The laminate is loaded to a strain e x in the x-direction. (A 
positive value of refers to tensile strain.) The associated strain in 

the y-direction is -v^ E x > where v lam is the base laminate Poisson's 
ratio. 
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The plan view of an assumed elliptic delamination is shown in figure 2(a). 
The set of intact laminae delaminated from the laminate is referred to as the 
"sublaminate ;** the remaining laminate is referred to as the base laminate. 

The principal axes of the sublaminate are x' and y f ; the corresponding 
semiaxes lengths are a and b. The angle between the axes x' and x is 
0, which is referred to herein as the sublaminate angle. The sublaminate is 
assumed to be made up of N laminae; a is the fiber angle of a lamina 
measured relative to the x-axis (see figure 2(a)). The sublaminate thickness 
is h and is assumed to be small compared to the base laminate thickness. 
Therefore, inplane displacements around the sublaminate boundary can be calcu- 
lated from the base laminate inplane deformations. Furthermore, the sublami- 
nate lateral dimensions (a and b) are assumed to be relatively large compared 
to h, and, hence, thin plate linear buckling theory is assumed to be valid. 

The buckled shape of a sublaminate is shown in figure 2(b). The trans- 
verse deflection w is measured from the sublaminate raid-plane. The trans- 
verse displacement and slopes are zero along the sublaminate boundary. The 
present analyses assume that the sublaminate buckles outward from the base 
laminate, as shown in figure 2(b). Furthermore, effects of higher modes and 
inward buckling of the sublaminate are neglected. The in-plane forces acting 
on the sublaminate due to the laminate strain were calculated from lami- 

nation theory [9] and are shown, schematically, in figure 2(c). Even though 
the base laminate is under uniaxial stress, the local stress distribution in 
the sublaminate is generally biaxial due to either a difference in the 
Poisson's ratio of the sublaminate and the base laminate, or to sublaminate 
material anisotropy. Note that the buckling initiation problem is reduced to 
that of an idealized elliptic sublaminate subjected to a set of inplane 
boundary forces. 
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ANALYSIS 


Two analyses were used to calculate the buckling strains of the sublami- 
nate: Finite-element (F-E) and Rayleigh-Ritz (R-R) . Both analyses were based 

on the adjacent equilibrium buckling criterion (see for example ref. 10). The 
finite-element analysis was performed using the STAGS [8] computer code. The 
F-E analysis, although versatile and accurate, is a relatively expensive 
computer analysis, especially for eigen value problems involving a large 
number of degrees of freedom. Therefore, an additional simplified analysis 
based on the Rayleigh-Ritz technique was developed. The R-R analysis is known 
to give accurate results for specially orthotropic sublaminates [11]. The 
term special orthotropy refers to an orthotropic sublaminate, one of whose 
principal axes is aligned with the applied load direction. In this study, the 
R-R method was used for both specially orthotropic and anisotropic sublami- 
nates. The accuracy of the analysis was assessed using the STAGS finite- 
element analysis. The two analyses are described in the following sections. 

The Finite Element Analysis 

The general purpose finite element program STAGS [8] was used in the 
present study. An eighteen degree-of-f reedom triangular plate element having 
three displacement and two rotational degrees-of-f reedom (d.o.f.) at the 
corner nodes and one rotational d.o.f. at each mid-side node was employed for 
the analysis. Specially orthotropic cases were analyzed using only a quarter 
of the sublaminate whereas generally anisotropic cases were analyzed using the 
full sublaminate idealization. The two idealizations are shown in figure 3. 

The following boundary conditions were used in the analyses: (1) Quarterplate, 
v = w, y , = 0 along Ok, u = w, x , = 0 along OB, and w = w, x , = w, y , = 0 
along AB; (2) Full plate, w = w, x , - w, y , = 0 along the entire outer bound- 
ary. The variables u, v, and w are the displacements in the x' ,y’ and 
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z directions, respectively, and the comma indicates partial differentiation 
(e.g., w, , = r)« For a specified base laminate axial strain e , the sub- 

laminate boundary forces (see fig. 2(c)) were input as initial stresses in the 
finite element analysis. 

Convergence was studied to determine the degree of mesh refinement 
required for specially orthotropic and generally anisotropic sublaminates. A 
quarter plate sublaminate of 8x12 (i.e., eight equal parts along the x'-axis 
and twelve equal parts along the circumference) mesh is shown in figure 3(a). 
Calculated buckling strains using 4x6, 8x12, and 16x24 meshes for a unidirec- 
tional 0 deg sublaminate indicate that the buckling strain for 8x12 mesh 
differed by less than 1% from that for 16x24 mesh. Thus, the 8x12 mesh was 
employed for all isotropic and specially orthotropic sublaminates. A similar 
convergence study for a unidirectional with a = 45 deg sublaminate (which is 
anisotropic) indicated that the 8x32 mesh (see fig. 3(b)) predicts buckling 
strain within 1 % from that of the more refined mesh. Thus, 8x32 mesh was 
employed for all anisotropic cases. 

The Raleigh-Ritz Analysis 

Even though the finite element method is versatile it is expensive, 
especially for eigenvalue analysis of problems with a large number of degrees 
of freedom. Therefore a simple Raleigh-Ritz analysis, based on the Trefftz 
criterion [10] , is presented to calculate buckling strains of elliptic sublam- 
inates. The procedure consists of three steps: (1) selection of a kinemati- 

cally admissible transverse displacement function; (2) calculation of the 
total potential energy; and then (3) application of the Trefftz criterion [10] 
to yield eigenvalue equations. The eigenvalue equations were solved numeri- 
cally for different sublaminate configurations and stacking sequences. The 
governing equations are derived below. 
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Strain energy of sublaminate.— The strain energy U of the buckled 
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where the D's are the flexural stiffness constants for an anisotropic sub- 
laminate [9]. The transverse deflection w of the deformed sublaminate is 
assumed to given by 


w = 1 1 


-w-fli 


C o + C 1 x 


• 2 + C 0 y 


,2 


( 2 ) 


where Cq, Cj, and C 2 are generalized displacements. The w function 
satisfies zero deflection and slope conditions along the boundary of the 
sublaminate. Substituting eqn. (2) in eqn. (1) and performing the necessary 
differentiation and integration, eqn. (1) becomes 


U = {C} T [K] {C> (3) 

where {C} T = {C Q , Cj , C 2 > and the stiffness matrix K is 

[K] = DhEKjl ] + D 22 [K 2 ] + 2[K 3 ](D 12 + 2D 66 ) (4) 


The matrices , K 2 and K 3 are defined as follows: 


8 



[K,l 

a 


4 0.5 

1.35 
Symmetric 


0.5 

0.05 

0.15 


[K ] = — 
2 J 

D 


4 0.5 

0.15 
Symmetric 


0.5 

0.05 

1.35 


HU = 


Gab 


8 1 

1.1 

Symmetric 


1 

0.5 

1.1 


Since the assumed w deflection function is symmetric, coefficients of Djg 
and £>26 vanish in the process of differentiation and integration; hence, 
these terms do not appear in eqn. (4). 

Potential energy of applied loads. - The potential energy V of the 
stress resultants is given by [9] 


V 2 f_ f, ( N X*Kc') 2 + N y'( W »y') 2 + 2N x «y'( W » x « W »y')j dX ' ^ 


where the inplane stress resultants N i, N v t, and N i i are obtained from 

A y x y 

lamination theory [9] as follows 


N X’ = A 11 e x' + ^2 V + e x’y * 

N y' = A 12 e x’ + A 22 e y' + A 26 e x'y' 

N x'y' = A 1G e x' + A 26 V + A 66 e x'y’ 

The sublaminate strains e , , e , and e . . 

x 1 ’ y' x y 

the base laminate strain e x , the Poisson's ratio 
angle 0 as follows 


( 6 ) 

are expressed in terms of 
V lam’ and t * le 8u blaniinate 
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( 7 ) 


2 2 

e , = (cos 0 - v, sin 0)e 
x *• lam J x 

2 2 

V " ( sln 9 - v latn “ s 9 ) C x 

Vy' - - (1 + " la J «in 26 = x 

Substituting the expression for w from eqn. (2) into eqn. (5) and performing 
the necessary differentiation and integration, eqn. (5) simplifies to 

V = {C} T [K g ]{C}e x (8) 

The geometric stiffness matrix K g is symmetric. The elements in the upper 
half of the [K g ] are defined as follows 

2 

r)[( A ll " v lam A 12 ^ + (§) ( A 12 " V lam A 22^] 

- [ a i 2 " v iam A 22 J 
\ t A ll “ v iam A 12 ^ 

(l)[ 3 ( A ll " V lam A 12 ) + (f) ( A 12 " V lam A 22^] 

2 

(b)[(l) ( A ll “ V lam A 12 ) + 3 ( A 12 “ V lam A 22^] (9) 

The shear-extension coupling terms A 16 and A 2 g do not appear in eqn. (9) 
because the w displacement function was symmetric. 
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The solution of eqn. (11) results in three eigenvalues which correspond to the 

three buckling modes. The lowest absolute eigenvalue corresponds to the first 

buckling mode and is referred to as the buckling strain e of the 

c 

sublaminate. 

To illustrate how the material anisotropy and the sublaminate size influ- 
ence the buckling strain, eqn. (11) is simplified to that for a one term 

solution. This corresponds to an assumed deflection function 

2 

C Q . Therefore, substituting Cj = C£ = 0 in 
eqns. (3) and (5) and applying the Treffetz criterion yields 



X cl 



( 12 ) 


where e is the one term buckling strain. Equation (12) shows that the 
X cl 

sign of e depends only on the sign of the denominator, which is a func- 
X cl 

tion of the size and the in-plane stiffness coefficients of the sublaminate, 

and the laminate Poisson’s ratio If the denominator is negative, the 

sign of the buckling strain e is positive. This indicates that the 

X cl 

sublaminate may buckle when the base laminate is loaded in tension. This 
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phenomenon may appear unusual, but it is caused by the mismatch of Poisson* s 
ratio between sublaminate and the base laminate. For example, if a quasi- 
isotropic base laminate having a Poisson's ratio of 0.3 is stretched to a 

strain e , it contracts laterally (in the y-direction) by -0.3 e . A delami- 
x x 

nated 90 deg fiber sublaminate, which has a very low Poisson's ratio (say 
0.03), contracts only by -0.03 The unequal contraction in the y-direction 

introduces a net compressive strain of -0.27 on the sublaminate. Such com- 

pressive strain acting in the fiber direction of the 90 deg fiber sublaminate 
may cause local buckling, even though the base laminate is loaded in tension. 

Equation (12) also shows that the sign of the denominator depends on the 
aspect ratio b/a of the sublaminate. The effect of sublaminate material 
properties and aspect ratio on buckling is discussed further in the next 
section. 

RESULTS AND DISCUSSION 

Buckling strains of a surface sublaminate in a quasi-isotropic graphite/ 
epoxy composite laminate were calculated using the finite-element (F-E) and 
the Rayleigh-Ritz (R-R) analyses. The R-R analysis was compared with F-E 
analysis for typical cases. Buckling strains for four-ply sublaminates were 
calculated for various sublarainate sizes, orientations, fiber angles and 
layups. The plot of buckling strains versus sublaminate properties is 
referred to as the buckling curve. Material properties used in the analyses 
are given in Table 1. 

Isotropic sublaminate. - Figure 4 compares buckling strains calculated 
from the F-E and R-R analyses for an isotropic (material properties for 
aluminum were used), elliptical sublaminate for different values of b. The 
F-E results are shown by symbols and the R-R results by a solid line. The two 
analyses agree very well. The compression buckling strain decreases in 


12 



magnitude raonotonically as the sublaminate width, b, increases and asymptoti- 
cally approaches the strip solution (b » a). 

Specially orthotropic sublaminates .- Figure 5 shows buckling curves for 
0 deg and 90 deg unidirectional sub laminates . The 0 deg unidirectional sub- 
laminate buckles only by compressive remote strain (as was the case for the 
isotropic sublaminate). The buckling strain asymptotically approaches a 
constant value with increased b. The 90 deg unidirectional sublaminate, as 
already explained with reference to eqn. (12), can buckle under tension and 
compression remote strain. The buckling caused by tension remote strain is 
referred to herein as tension buckling and buckling due to compression remote 
strain is referred to as compression buckling. Results in figure 5 and an 
examination of buckling modes (not reported here) Indicate that the first mode 
of buckling for a 90 deg unidirectional sublaminate only occurred for remote 
tension for b < 25 mm and only for remote compression for b > 64 mm. Both 
compression and tension buckling of a sublaminate are possible for values of 
b between 25 and 64 mm. Higher modes of buckling involve both positive and 
negative lateral displacements. The base laminate would interfere with these 
negative lateral displacements of the sublaminate. To account for this inter- 
ference requires the solution of a nonlinear contact problem, which is beyond 
the scope of the present paper. Hence, higher mode buckling strains are not 
reported here. 

The tension buckling strain for b - 12.5 mm is 0.0051, which is higher 
than the splitting threshold for unidirectional graphite/epoxy laminate 
(transverse tensile strength - 40 MPa, ref. 12). Therefore, the 90 deg 
unidirectional sublaminate for b < 19 mm may split before buckling. The 
splitting would reduce the buckling strain. The two analyses agree very well 
with each other for both sublaminates. 


Anisotropic sublaminates. - The lower half of figure 6 shows the effect of 
fiber angle a on buckling strain for unidirectional fiber, circular sublami- 
nates. The compression buckling strain increases with increasing fiber angle. 
But at fiber angles of 75 deg < a < 90 deg, the sublaminate buckles also 
under tensile laminate strain. The F-E and R-R analyses agree well for 

0 < a < 55 deg and 80 deg < a < 90 deg. For fiber angles 55 deg < a < 80 deg, 

the R-R analysis overestimates the buckling strain, because the bending 
twisting (D 16 and D 26 ) and shear-extension (A 16 and A 26 ) coupling terms 
were not included. However, for other fiber angles the influence of these 
terms is less significant and hence the R-R results agree with the F-E results. 

The upper portion of figure 6 also shows the effect of fiber angle on 

major Poisson's ratio (v . ,) of the sublaminate. The Poisson’s ratio of the 

x y 

base laminate is also shown in figure 6 for reference. Comparison of the 
buckling curve and Poisson's ratio curves indicate (1) the buckling strain of 

the sublaminate is compressive and is relatively small for sublaminate 

Poisson's ratio greater than the base laminate; (2) buckling strain increases 
in magnitude as the sublaminate Poisson’s ratio decreases and falls below the 
laminate value; and (3) the first mode of buckling changes from compression to 
tensile buckling for very low sublaminate Poisson's ratio. 

Figure 7 shows the effect of fiber angle on the buckling strain for three 
unidirectional elliptical sublaminates with different aspect ratios. The 
length "a" was held constant (a - 25.4 mm); the widths b » a/2, a and 2a. 
Results for the circular sublaminate in figure 6 are repeated here. Buckling 
strains for b/a =0.5 were compressive for a < 60 deg and became tensile 
for a close to 90 deg. In contrast, buckling strains of b = 2a sublami- 
nate remained compressive for all fiber angles. Buckling strain magnitude 
decreased monotonically with increase in sublaminate width, b. 
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Oblique sublaminate. - The effect of fiber angle on the buckling of oblique 
sublaminates is shown in figure 8 . Buckling curves are shown for unidirec- 
tional sublaminate angles 0=0, 15, 30, 45, 60, and 90 degrees. The sublami- 
nate size is a = 25.4, b = 50.8, and h = 0.51 mm. Figure 8 shows that for 
fiber angles 0 < a < 38 deg, buckling strain decreases slightly with increasing 
sublaminate angle. The trend is reversed for a > 38 deg. For 9=0 and 
15 deg the buckling strain is compressive for all fiber angles; whereas for 
30 deg < 0 < 90 deg the buckling strain changes to tensile for a > 70 deg. 

Figure 9 shows the effect of sublaminate angle on buckling strain for 
[°l 4 > [90]^, [0/90] g , [±45] g , and [0/±45/90] sublaminates. The figure shows 
that the buckling strains for [0/90] g , [±45] g , and [0/±45/90] sublaminates are 
generally bounded by those of the [0]^ and [90]^ sublaminates. The buckling 
strain for a [ 0 ]^ sublaminate decreases with increasing sublaminate angle 
0; for a [90]^ sublaminate the buckling strain increases. Furthermore, 
buckling strains of [0 ] 4 , [0/90] g , [±45] s , and [0/±45/90] sublaminates are 
quite small (less than 0.0015) and compressive for all sublaminate angles. 

Criterion for delamination growth direction .- For the configuration being 
studied, there are no interlaminar stresses around the delamination boundary 
until the delaminated region buckles. Whether the delamination grows after 
buckling depends on many parameters, including sublaminate size, sublaminate 
in-plane and flexural stiffness, and the interlaminar fracture toughness of 
the material. If we assume the delamination will grow, the question remains: 
"What direction will the delamination grow?" The possibility that the direc- 
tion of initial growth for a circular delamination could be predicted based on 
changes in the buckling strain was investigated. The delamination was assumed 
to be Initially circular (a = b = 25.4 mm). Then, it was assumed to grow in 
the load (x) direction by 50% (a = 38.1, and b = 25.4 mm) or transverse to the 



load direction by 50% (a = 25.4, and b = 38.1 mm). Buckling strains for the 
three configurations were calculated for [O]^, [0/90] s , [±45] g , and [0/±45/90] 
sublaminates. Figure 10 shows the calculated results and the experimentally 
observed growth directions [1,2,13]. The [0]^ sublaminate buckling strain is 
lowest when it is elongated in the load direction. In tests reported in 
reference 13 the delamination also grew in the load direction, accompanied 
with splitting. On the other hand the buckling strains for [0/90] g , [±45] g , 
and [0/±45/90] sublaminates are lower when the delamination is elongated 
perpendicular to the load direction, which is also the experimentally observed 
growth direction [1,2]. Apparently, the direction of elongation of the sub- 
laminate that has the lowest buckling strain will be the initial direction of 
growth of a circular delamination. More tests and analyses are needed to 
verify this interpretation. 

CONCLUSIONS 

Finite-element (F-E) and Rayleigh-Ritz (R-R) analyses were used to 
predict the buckling strain of a localized delaminated region in a quasi- 
isotropic laminate. The R-R analysis, developed using the Trefftz criterion, 
is simple and inexpensive compared to the F-E analysis. The accuracy of the 
R-R analysis was evaluated by comparing with the F-E analysis. Calculated 
buckling strains from the R-R analysis agreed with the F-E results except for 
highly anisotropic sublaminates; the neglect of the bending-twisting (D^g and 
D 2 g) and the shear-extension (A^g and A 2 g) coupling terms caused significant 
errors for these sublaminates. However, the region where the R-R method is 
invalid was small. A parametric study was performed for a four ply sublami- 
nate to assess the influence of material anisotropy, and sublaminate angle 
shape and layup on buckling strains. The laminate was assumed to be made up 
of graphite/epoxy composite. This study led to the following conclusions: 
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(1) Even though the base laminate is in a uniaxial stress state, the sub- 
laminate is generally in a biaxial stress state either due to a mismatch 
of Poisson 1 s ratio or to material anisotropy. 

(2) The biaxial stress states in some unidirectional composite sublarainates 
cause tensile buckling when the sublaminate is elongated in the load 
direction and fibers are nearly perpendicular to load direction. 

(3) Compressive buckling strain for unidirectional sublarainates increases with 
the increase in the angle between the load and the fiber directions. 

(4) The buckling strains for [0/90] g , [±45] s , and [0/±45/90] sublarainates are 
generally bounded by those for [0]^ and [90]^ sublaminates. 

(5) The direction of elongation of the sublaminate that has the lowest buckling 
strain correlates with the initial growth direction of a circular 
delaraination. 
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Table 1.- Material Elastic Properties 


Material 

Modulus, GPa (MSI) 

Poisson ratio 

E * 

E t 

G *t 

V At 

Aluminum 

68.95 

(10.0) 

68.95 

(10.0) 

26.32 

(3.82) 

0.31 

Graphite/epoxy 

131.0 

(19.0) 

13.0 

(1.89) 

6.4 

(0.93) 

0.34 
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Figure 1.- Shapes of delaminations. 
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(c) FREE BODY DIAGRAM 
OF THE SUBLAMINATE 


Figure 2.- Nomenclature for a buckled sublaminate. 










Figure 4.- Aluminum sublaminate: effect of sublaminate width on buckling strain, 

(a = 25.4 mm, h - 0.51 mm) 
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Figure 5.- Specially orthotropic sublaminate: effect of sublaminate width (b) on 

buckling strain. (a = 25.4 mm, h = 0.51 mm, graphite/epoxy) 
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Figure 6.- Circular sublaminate: effect of fiber angle on buckling strain, 

(a = 25.4 mm, h « 0.51 mm, graphite/epoxy) 
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Figure 7.- Elliptical sublaminate: effect of fiber angle on buckling strain, 

(a = 25.4 mm, h = 0.51 mm, graphite/epoxy) 
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Figure 8.- Unidirectional oblique sublaminate: effect of fiber angle and sublaminate 

angle (0) on buckling strain. (a = 25.4 mm, b = 50.8 mm, h = 0.51 mm, 
graphite/epoxy) 
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Figure 10.- Buckling strains for three sublaminate configurations. 

(h = 0.51 mm, graphite /epoxy) 
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